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Abstract 

We consider the time-independent Green function for the Schrodinger operator 
with a one-particle potential, defined in a d-dimensional domain. Recently, in one 
dimension (ID), the Green's function problem was solved explicitly in inverse form, 
with diagonal elements of Green's function as prescribed variables. In this paper, the 
ID inverse solution is used to derive leading behavior of Green's function close to the 
domain boundary. The emphasis is put onto "universal" expansion terms which are 
dominated by the boundary and do not depend on the particular shape of the applied 
potential. The inverse formalism is extended to higher dimensions, especially to 
3D, and subsequently the boundary form of Green's function is predicted for an 
arbitrarily shaped domain boundary. 
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1 Introduction 



For a connected finite domain, the eigenvalue spectrum of an operator with a boundary 
condition depends strongly on the shape of the smooth boundary. A rigorous analysis 
of the eigenvalue density was first done by Weyl 1 for the 3D Laplacian in the asymp- 
totic limit of an infinite volume (at short wavelengths of eigenstates). The next surface 
and curvature contributions to the eigenvalue density were derived for the Dirichlet and 
Neumann boundary conditions first in 3D, 2 and then, by using the path integral method, 
in 2D. 3 In this connection, Kac posed a question whether there is a one-to-one corre- 
spondence between the boundary shape and the corresponding eigenvalue spectrum of 
the Laplacian. The answer is negative: 4 there exist nonisometric pairs of 2D flat shapes 
which are isospectral. 5 ' 6 An important progress in calculations was made within a time- 
independent Green function method in Ref. 7. Here, a multiple reflection expansion was 
used to get Green's function for an arbitrary domain, in the limit of small wavelengths. 
This technique, applied to physical problems like electromagnetic field in a cavity, 8 evo- 
qued numerous research activities 9 . Its extension to the Schrodinger operator 10 resulted 
in semiclassical expansions for quantum mechanical system valid for finite domains 11 ' 12 . 
The influence of spatial confinement on the energy spectra of simple quantum systems, 
like harmonic oscillator or hydrogen atom, was investigated in many other works (see e.g. 
Refs. 13-17). 

This paper deals with local effects of a domain boundary on Green's function of a 
Schrodinger operator. The topic is studied within an "inverse" formulation for Green's 
functions. The motivation for this formulation comes from the seemingly unrelated 
density-functional theory. The latter is based on the famous Hohenberg-Kohn unique- 
ness theorem, 18 having a practical realization in the Kohn-Sham separation ansatz. 19 In 
the density-functional theory, the usual "direct" formulation, i.e. the ground-state ex- 
pectation value of an observable in terms of some prescribed external fields, is replaced 
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by the "inverse" one, with the particle density as the controlling variable. In a series of 
works, 20-22 we have developed an infinite gradient series expansion of the kinetic-energy 
density functional for ID noninteracting Fermion systems in the ground state by applying 
the Green's function method. By definition, the expansion in potential gradients and 
correspondingly in density gradients does not allow us to consider discontinuities, say at 
hard particles or extended regions. The "inverse" formulation of Green's functions arised 
as a natural mean for their explicit representation in any situation. 

To introduce the inverse format for Green's functions, let us first specify notation. We 
consider a d-dimensional domain Q of points r = (x±, x 2 , ■ ■ ■ , xj), infinite (r e R d ) or finite, 
bounded by hard walls at the domain boundary dVL. Selecting units so that h 2 /(2m) = 1 
with m being the particle mass, the general one-particle Hamiltonian formulated within 
the domain Q is given by 

H = — A r + u(r). (1.1) 

Here, u(r) is an arbitrary external potential, continuous inside Q, which does not include 
the infinite potential due to the presence of domain walls at dQ. The time-independent 
Schrodinger equation reads 

H^(r) = A fc ^ fc (r) (1.2) 

with the Dirichlet b.c. ipk{ Y ) — at the domain boundary r G dQ. ipk{ r ) an d A& are 
respectively the eigenstates and the eigenvalues of Hamiltonian H, and index k can be 
either discrete (localized eigenstates) or continuous (extended eigenstates). The time- 
independent Green's function is defined by 23 G z = (z — H) _1 , where z is a complex 
variable with components A = Re(z) and s = lm(z). G z is an analytic function of z 
except of those points on the real z-axis which correspond to the eigenvalues of H: it 
exhibits simple poles at discrete eigenvalues of H and a branch cut along fragments of 
the real z-axis which correspond to the continuous spectrum of H. In the latter case, the 
discontinuity of the transverse limits = lim s ^ + Ga±is yields the density of states at 
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A. In the vector-space representation, 

G , (ri0= £*M5fcOO (L3) 

k Z — Ak 

Green's function satisfies the two-point differential equation 

[A r + z-u(r)]G z (v,r') = 5(r-r') (1.4) 

with the b.c. at the domain boundary, G z (r, r') = at r e cftl We notice that in one 
dimension (ID), diagonal elements of Green's function, denoted by 



n,(x) = E TT ' (L5) 



are well defined finite quantities. This is no longer true in dimensions > 2 where the 
equal-argument Green's function diverges. One can avoid this divergence by redefining 
new one-point quantities, which choice is crucial from the point of view of this work. In 
3D, which is of practical interest, the finite quantity 

u z (t) = Iim47r|-G' 2 (r,r / ), d = 3, (1.6) 

will be especially useful. The function v z remains informative since the differentiation of 
Green's function with respect to z keeps the poles of the latter. One can reconstruct the 
original Green's function by using the relation 

f ds^(r)=47r(r|— ^ + i|r>, d = 3. (1.7) 
Jo z - H H 

The ^-independent term in (1.7), which subtracts the diverging part of the diagonal 
Green's function element, has no effect on the density- functional formalism developed in 
Refs. 20-22. 

The above format is the direct one: for a given constraining domain Q and a given 
external field w(r), find the corresponding Green's function. For continuum ID space 20 ' 21 
and for simply connected lattice structures, 22 the direct problem was replaced by the 



inverse one with diagonal elements of Green's function as controlling variables: (i) the 
external field in terms of the diagonal elements of Green's function (the inverse profile 
relation); (ii) off-diagonal Green's function elements in terms of the diagonal ones. The 
ID inverse profile relation was found explicitly in the form 
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z - u(x) = — + - 
An 2 z (x) 4 



n 'z( X ) 

n z (x) 



1 n" z {x) 
2n z (xY 



with the standard b.c. n z (x) = at x G dVl. A prime means the derivative with respect 
to the argument. Eq. (1.8) with b.c. determines n z (x) uniquely up to the sign. When 
compared to the basic equation (1.4) of the direct format, it has the virtue of being a 
one-point differential equation without the appearance of a 5-function. As concerns the 
off-diagonal elements of Green's function, they can be expressed in terms of the diagonal 
ones as follows 



G z (x, x') = [^(rr)] 1 / 2 [n^rr')] 1 ^ 2 exp sign(x — x ) 

Jx' 



r x ds 



(1.9) 



2n z (s)_ 

As was already mentioned, in this paper we use the inverse formalism to study behavior 
of the Green's function close to the confining domain boundary. In the leading and some 
higher orders of the distance from the boundary, the Green's function is shown to exhibit 
" universal" terms which are dominated by the boundary and do not depend on the applied 
potential. These boundary terms are derived first in ID, and subsequently predicted in 
3D for an arbitrarily shaped domain boundary. The results are exact at all wavelengths. 

The paper is organized as follows. In Sec. 2, we derive in the inverse format the 
boundary form of ID Green's function. Before extending the inverse formalism to higher 
dimensions, in particular to 3D, we first treat, in the direct format with the zero potential, 
two special cases which can be treated exactly as one-dimensional. Systems that are 
stratified to only one Cartesian coordinate (Sec. 3) justify the choice of the one-point 
quantity v z (1.6). Spherically symmetric systems (Sec. 4) reveal unexpected logarithmic 
terms in the boundary expansion of Green's function. In the latter by-product 



of the formalism, we derive general behavior of Green functions close to the center of 
radial symmetry. A technique for systematic construction of the inverse profile relation 
in an arbitrary dimension and for any potential is developed in Sec. 5. As a consequence, 
the (universal) boundary terms in 3D Green's function are suggested for an arbitrarily 
shaped domain boundary. 



2 ID Green's function 



We start with ID Green's function G z (x, x') defined in a domain f2 = (0, X), and look for 
its behavior close to a wall, say as x, x' — > 0. The potential u(x) is assumed to be regular 
at x — 0, and its Taylor expansion for x > is written as 

1 & n u(x) 



oo 



U ( X ) = Yl U n% n , 
n=0 



U n = 



x=0 



(2.1) 



n\ dx n 

The first coefficient u only shifts z — > z — u , so without any loss of generality it is set 
equal to 0. 

Let us suppose that the diagonal element of Green's function n z (x), satisfying the b.c. 
n z (0) = 0, also has the Taylor series expansion around x — 0, 

1 & n n z (x) 



n 



(x) ^ ' c n x , c r , 



n=l 



x=0 



(2.2) 



n\ dx n 

and insert this expansion into the exact inverse profile relation (1.8). The requirement of 
the vanishing of the prefactors attached to °, etc., implies the following sequence 

of equations for the coefficients {c n }: 



(ci) 2 -l = 0, 



c 2 



[ci? - 1 



0, 



1zc\ + 3C3 = 0, 



(2.3a) 
(2.3b) 
(2.3c) 



etc. The first Eq. (2.3a) tells us that K(0)] 2 = I. Actually, we observe in ID examples 



with boundaries, like free particle or harmonic oscillator in bounded well, that 

n'M = -1. (2.4) 

This means that n z goes to zero at the boundary from below. The universal slope —1 
depends neither on details of the applied potential nor on the b.c. at the opposite bound- 
ary X. Relation (2.4) also tells us that the functional series, which determines n z (x) via 
formula (1.5), is not uniformly convergent. In the opposite case, we could differentiate 
this functional series term by term with respect to x, which leads to the contradiction 
n' z (0) = 0. The second Eq. (2.3b) is fulfilled identically, leaving the coefficient c 2 unspec- 
ified: the value of Ci{z), which depends on z and {«„}, is fixed by the requirement of 
the vanishing of n z [x) at the opposite X-boundary. The third Eq. (2.3c) determines the 
coefficient c 3 , which, being the function of only z (trivially shifted by u ), is universal in 
the same sense as c\. The next coefficients c 4 , c 5 , • • • are analytic functions of both {u n } 
and C2, which confirms the adequacy of the analyticity assumption (2.2). We conclude 
that 

n z {x) ~ —x + ci[z)x 2 + -^-x 3 asi-> + . (2.5) 

When the two points x and x' are close to the boundary at 0, the leading terms of 
G z (x,x') can be evaluated by inserting (2.5) into (1.9), with the result 

G z (x, x) ~ — x < + c 2 {z)xx' + -x < (x\ + 3rr?^ as both x, x — > + . (2.6) 

Here, we have introduced the standard notation x< = min{x, x'} and x > = max{i,a;'}. 

3 d-dimensional Green function with a ID confine- 
ment to potential 

There exists a large family of rf-dimensional models whose Green's function can be ex- 
plicitly expressed in terms of a related ID Green's function. The family is defined by 
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potentials that are stratified in one dimension, i.e. depend on only one Cartesian coor- 
dinate of r = (xi, ■ ■ ■ ,Xd), say X\ = x G Q x , u(r) = u(x). For simplicity the subspace 
of vectors r_i_ = (x 2 , ■ ■ • ,Xd) perpendicular to the rr-axis will be infinite, Q = Q x <g> i2 d_1 . 
Using the separation-of-variables treatment with plane waves kj_ = (k 2 , ■ ■ ■ ,kd) in the 
perpendicular subspace, the <i-dimensional Green's function is expressible as 

where G 1D is ID Green's function associated with the potential u(x) and the domain fl x . 
In d — 3 dimensions, passing to polar coordinates, formula (3.1) takes form 

1 r°° 

Gf (r,r') = — / dA; kG™ k2 (x,x')J (kp ± ). (3.2) 
2ir Jo 

Here, p± = \r± — r' ± \ is the perpendicular distance of points r, r' and J denotes the Bessel 
function of the first kind. 24 

Let the 3D domain f2 be the half-space (0, oo) (g> R 2 . The corresponding ID Green's 
function G^_ k2 (x, x') in (3.2), formulated within fl x = (0, oo), can be expanded in x and x' 
near the rectilinear hard wall at using the previously derived formula (2.6). It might be 
tempting to insert this expansion directly into (3.2), but already the leading term ~ —x < 
implies a diverging integral over k. This indicates that a renormalization procedure is 
needed for ID Green's function. 

We first assume that p± ^ 0, i.e. the two points r and r' do not lie on the same line 
perpendicular to the surface of the wall. G 1D can be split into two parts 

+ 5G™(x,x'), x,x'>0, (3.3) 

where one assumes that lm(y/z) > in order to ensure the regularity at \x — x'\ — > oo. 
The first part, which is nothing but ID Green's function of free particle in the half-space 
x > 0, contains all universal terms of the small-x, x' expansion (2.6). The remainder does 
not contain the universal terms and behaves like 

xx 1 + 0[xx\x 2 + x'% (3.4) 
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G™(x,x>) 



5G™(x,x')= cf{z) + i^-z 



where the model-dependent c^(z) is defined by (2.5) or (2.6). Substituting (3.3) and 
(3.4) in the basic formula (3.2), one gets 



where 



47T 



exp(V£|r-r'|) exp [iy/zyj(x + x') 2 + p z ± 



r — r 



[x + x') 2 + p\ 



+ SG™ (r,r% (3.5a) 



^ D ( r - r ') = -^f dSS 



C 2 D (^s 2 ) + iv^s] h (V^P±Vs 2 - l) xx' + O [xx'(x 2 + x' 2 ) 

(3.5b) 

It can be easily shown that clp(z) + iy/z ~ 0(1/ z) as \z\ — > oo, so the integral in (3.5b) 
has the necessary convergence property. Consequently, the total G^ D (x, x'\ p± ^ 0) oc xx' 
does not exhibit universal expansion terms near the boundary, as is intuitively anticipated. 

When p± = 0, i.e. the two points r and r' lie on the same line perpendicular to the 
surface of the wall, relation (3.2) can be transformed to 



Gf (x,x'; P± =0) =Gl D (x,x';p ± = 0) + ^- [* ds Gf(x,x'). 

Air Jo 



(3.6) 



It is straightforward to show that the small-x, x' expansion of Gq D is determined by the 
free-particle limit plus a model-dependent term ~ 0(xx'). From (2.6), Gl B (x,x') = 
—x K + 0(xx'). Substituting this in (3.6), one observes that 



GT(x,x'-p ± =0) = 



x < 



-t—x< + O(xx'). 



(3.7) 



2tt\x — x'\ (x + x') 4ir 

The first two leading terms on the rhs of (3.7) are universal, the applied potential con- 
tributes starting from the term of order 0(xx'). 

We add that the 3D one-point quantity u z , introduced in (1.6), equals the diagonal 
element of ID Green's function associated with the potential u(x), and therefore satisfies 
the inverse relation of type (1.8), 
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Z - u(x) = -- z + - 

4:U 2 (x) 4 
and exhibits the small-x expansion (2.5). 



v' z (x) 



VAX) 



1 <(*) 

2 v z (x) ' 



(3.8) 
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4 Radial Green's function 



We now consider a 3D quantum system with radial symmetry, confined to the domain 
^ = {|r| < R} (radius R may be either finite or infinite). The external potential u(r) 
depends only on the magnitude r of r. The radial problem can be reduced to ID by 
factoring out the angular dependence of the Hamiltonian eigenfunctions in terms of the 
spherical harmonics. Let us introduce for each angular momentum quantum number 
I — 0, 1, 2, . . . the ID Green function associated with the Hamiltonian 



H(0 



d 2 1(1 + 1) 

+ 9 ' +u(r), 



dr 



(4.1) 



= (z — fiW) with zero b.c. at the origin r = and at the boundary r = R. The 
total 3D radial Green's function is then expressible as 

(2/ + i; 



Gl D lr 



1=0 



(4.2) 



where uj is the angle between r, r' and Pi denotes the Legendre polynomial of degree Z. 24 
The radial analogue of the ID inverse profile Eq. (1.8) for the diagonal elements of 



n 



(r) = G^(r,r), reads 



z ~ M ( r ) 3— 



4ni°(7 



1 

\2 4 



(r) 



lr£V) 
2 r$(r) ' 



(4.3) 



with the obvious b.c. n®(0) = n^\R) = 0. According to (1.9), the off-diagonal elements 
of G^ are expressible in terms of the diagonal ones as follows 



Gj)( r y)= [nj>(r)l 1/a \4\r') 



1/2 



exp 



sign(r — r') 



ds 



r' 2n i z\s). 



(4.4) 



4.1 Behavior close to the center 

Although the additional potential term 1(1 + l)/r 2 (I ^ 0) is singular at r = 0, the r — > 
analysis of Eq. (4.3) is similar to that of the ID inverse profile relation (1.8) close to a 



10 



boundary (see Sec. 2). Let us assume that the potential u(r) is regular at the origin. 
Inserting into (4.3) the Taylor expansion of around r = 0, 



oo 



n=l 



r=0 



(4.5) 



one gets a sequence of equations for the coefficients {c®} which implies 



4° M + $ o4 0) (^)r 2 - — — - - — -r 3 as r -> 0. (4.6) 

zKJ 21 + 1 l ^ 2KJ (21- l)(2/ + l)(2/ + 3) V ; 

The leading terms of the expansion of Green's function G®(r,r') in r and r' can be 

evaluated by inserting (4.6) into (4.4). The total 3D Green's function (4.2) can be obtained 

by applying the generating formula for Legendre polynomials, with the result 

Gf(v,Y')~-—± + — |r-r'| asbothr,r'^0. (4.7) 

2 47r|r - r'| Att 8tt 

Here, the potential u(r) and the presence of the boundary at r = R are reflected only in 

the model-dependent coefficient ^(z). 

4.2 Behavior close to the boundary 

We introduce a new variable, the distance from the wall x = R — r, and redefine u(r) — > 
u(x), n®(r) -> n(' } (x), Gf'(r,r') -> Gf'^i')- The total 3D Green's function (4.2) is 
expressible simply as 

Gf{x, x';u) = _ _ ^ g ^P,(oos,) G l"(,, (4.8) 

The "angular momentum" potential 1(1 + 1)/(R — x) 2 is an analytic function of x close 
to the boundary. Provided that u(x) is also regular at x — 0, it holds G®(x,x') 
as x, x' — > 0. Inserting this into (4.8), the 3D Green's function becomes proportional to 
the non-converging series X^ (2/ + 1)P^(cosc<j) and we face the same problem as in the 
previous case of the 3D rectilinear boundary (Sec. 3). 

Since we know how to solve the divergence problem for the rectilinear hard wall, our 
strategy is first to reproduce the result (3.2), valid for the 3D rectilinear wall, as the 
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R — > oo limit of the relation (4.8), and then to get the leading 1/R correction due to 
the curvature of the wall surface. We have treated the angular momentum potential 
1(1 + l)/(R — x) 2 as a perturbative series in 1/R, then used the standard Green's function 
perturbation theory for and finally performed a convenient continualization of the 
sum on the rhs of (4.8) setting I — kR — 1/2 (k fixed) and 

Piictxu) = Jo(k P± ) - kp±i ^ X ' ] Mk P± ) + O (^) . 

Here, p±, defined by p\ = |r — r'| 2 — (x — a;') 2 , is the analogue of the perpendicular distance 
|rj_ — r'jj in the case of the rectilinear wall. The final result is 

Gf (x,x';p ± ) = —J o dk kG™ k .{x,x')Uk Pi _) + -g z ( x , x >-p ± ) + o(—) , (4.9) 



where 



G z {x,x'; PA _) = ^±^1 T dkkGf_ k ,{x,x')Ukp A 

Z7T JO 



p±(x + x') f°° 



47T 

1 r°° 

+ 



/•OO 

/ dk fc 2 G^ fc2 (x,x')Ji(A;p ± ) (4.10) 
Jo 



- dkk 3 J (kp ± ) dyGl%(x,y)yG 1 ^(y,x'). 

7T Jo JO 



G 1D is the ID Green's function corresponding to the given potential u(x), with the zero 
b.c. at x = and the regularity condition at x — > oo, 

The first term in (4.9) is nothing but the exact result (3.2) valid for the R — > oo 
rectilinear hard wall. Its expansion for small x, x' coordinates was discussed in Sec. 3. 
The next term Q z {x, x'; p±)/R represents the leading correction due to the curvature of 
the wall surface. In the particular case of the zero potential, when G 1D is given by (3.3) 
with 5G 1D = 0, the small-x, x 1 expansion of Q z reads: 



xx'z 



Q z (x,x';p ± ^0) = —K {-iyG P± ) + 0\xx\x + x% (4.11) 

where K is a modified Bessel function; 24 

G z {x, x'- p ± = 0) = - XX ' - ^ In \-i^(x + x')} + O(xx'). (4.12) 
Atx[x + x'Y Ait 1 J 
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Worked-out examples with nonzero potentials tell us that for p± ^ the leading term 
oc xx' on the rhs of (4.11) depends on the applied potential, while there is evidence that 
for p± = the first two leading terms on the rhs of (4.12) are universal (i.e. independent of 
the applied potential) and the potential enters into the next term of order 0(xx'). From 
(4.12), the one-point quantity of interest v z has the following short-distance expansion 
from the wall 



where x = \/ y/z is the length parameter. We see that the nonzero curvature of the 
sphere surface induces a non-analyticity, namely divergence of the second and higher- 
order derivatives of v z with respect to x at the sphere surface x — 0. In the following 
section, we will prove the universality of the logarithmic term in the expansion (4.13) and 
derive its general form for an arbitrarily shaped domain boundary. 

5 d-dimensional Green's function 

Let us now consider a rf-dimensional quantum system with the general Hamiltonian (1.1) 
and Green's function satisfying Eq. (1.5). Representing the 5-function by 5(r — r') = 
/ exp [ik • (r — r')] dk/(27r) d , k = (k\, . . . , kd), we rewrite (1.5) as follows 




(4.13) 




(5.1) 



It is easy to derive the commutation rule 



\z - u(r) + Vr] _1 exp(ik 



exp(ik -r)[z- u(r) + (V r + ik) 2 ] 



i 



(5.2) 



Consequently, 



G z (r,r>) = J 



dk 



exp [ik 



r')] [z-u^ + iVr + ik) 2 ]- 1 



(2%) d 



I 



dk 



exp [ik 




z — u(r) — k 2 \ z — u(r) — k 2 



(5.3) 
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The standard expansion of the inverse operator in (5.3) then results in an infinite series 
1 1 r (V 2 + 2ik-V) + 



z — u(r) — k 2 z — u(r) — k 2 " z — u(r) — k 

where we have dropped the subscript r from V. Writing the series (5.4) formally as the 
sum J2 n «ri( r , — u ( r ) — k 2 ] n , we have at once 

~ a w (r,k) _ 1 1 (V 2 . 2jk , ^ a n (r, k) 

^[z - u(v) - k 2 ] n z-u(r)-k 2 z-u(v)-k 2K ' ^[z - u(r) - k 2 ] n ' 

(5.5) 

This differential equation implies the following recursion for the coefficients {a(r, k)}: 

«n = S n ,i ~ [V 2 a n _i + 2ik • Va n _i] 

-(n - 2) [2ia n _ 2 k • Vw + 2Va n _ 2 • Vw + a n _ 2 V 2 w] (5.6) 
—(n — 2)(n — 3)a n „3| Vw| 2 . 

Green's function (5.3) is expressible in terms of the coefficients {«„} as follows 

G ^ r '» = 1 ex » i* ' < r - r ')l | [.-°t)'-W - (5 ' 7) 

Introducing spherical coordinates in k-space, the integrals on the rhs of (5.7) can be 
evaluated by using the residue theorem in the complex A;-plane. 

It is straightforward to retrieve the ID profile relation (1.8) by using the above scheme. 
The calculations are more tedious in 3D. Generating the coefficients {«„} from the recur- 
sion (5.6), formula (5.7) gives, in the direct format, 

1 1 1 Au(r) 



liMr 



The inversion of (5.8) starting from the local-homogeneity reference z—u(r) = — l/[4u 2 (r)} 
results in the inverse profile equation 

z - u(r) = -—^ \l - |W 2 | 2 + 2u z Au z ] - £ ip n (r). (5.9) 



4i£(r) 



n=l 
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Here, the terms which involves just 2n z/ z -functions are grouped into y? n . In particular, 



<Pi 



(5.10a) 



i,j,k 



¥2 = H d A u z) 2 ( d kV z ) 2 + ^(d i d k u z )(d j d k u z )(d i u z )(d j u z ) 



A(dlv z )(d 2 i d j v z )(d j v z ) - 2{d 2 d j v z )(d j d k v z )(d k v z ) 



i,j,k 



-2(d i d j d k u z )(d i d j u z )(d k u z ) - 2(d i d j u z )(d j d k u z )(d i d k u z ) 



+{d i d j v z )\div z ) + (dfu z )(du z )(diu z ) 



i,j,k 



+ {d 2 i d k v z ){d 2 j d k v z ) - (didjd k iy z y 



(5.10b) 



etc. Note that the convergence of the direct series (5.8) is restricted to the high-energy 
part of the spectrum, z — u(r) >> 1, and to slowly varying potentials u(r). These are the 
attributes of extended states, and, indeed, v z given by (5.8) exhibits a branch cut along 
real z-axis. On the other hand, the inverse profile relation (5.9), being supplemented by 
the appropriate b.c. for u Z) holds in the whole complex z-plane, including simple poles 
induced by discrete localized states of the Hamiltonian spectrum. This is an important 
feature of the inverse format. 

As a check, we know from Sec. 3 that (5.9) must reduce to the ID profile relation 
(3.8) when both u(r) and v z {y) are stratified in one dimension. Under these circumstances 
<f 1 — ip 2 — ■ • • — 0, and indeed (5.9) becomes identical to (3.8). 

Like in the ID case, the 3D inverse relation (5.9) is very appropriate to describe the 
behavior of i^(r) close to the boundary dfl, where v z — > 0. Provided that u(r) is regular 
at r G dQ, the lhs of (5.9) is regular at r G d£l and so must be also the corresponding 
rhs. The cancellation of the leading diverging terms on the rhs of Eq. (5.9) is determined 
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exclusively by the terms oc 1/z/ 2 , i.e. by 

z - u(r) ~ --^rr [l " |V^| 2 + 2u z Au z ] . (5.11) 

This relation implies all universal terms of the //^-expansion around the boundary. 

Let us analyze relation (5.11) for the general potential u(r) and the general 3D domain 
fl with a smooth 2D boundary dfl of points r defined implicitly as follows 

dfl: 0(r o ) = Co. (5.12) 

The function is such that C has the dimension of length [for example, 0(r o ) = 
\Jxq + Uq + Zq(= R) for the sphere]. The set of points {r } G dfl can be parametrized by 
two curvilinear coordinates 9 and ip, r = r (9,ip), which are chosen to form an orthog- 
onal coordinate system. To every point r G fl, we attach the number C = 0(r) (with 
the dimension of length) and introduce the coordinate £ = C — C . The point r e dfl 
adjoint to a given r results as the intersection of the surface dfl with a curve which passes 
through r and simultaneously is perpendicular at every point r' to the surface 0(r') = C. 
The relationship between r and r reads 

oo 

r = r + £ u n (r )C, (5.13) 



71=1 



where 



= 2|V0(r o )l 4 I^d>(ro)dj0(r o ) 

(),<Art)) ^9 J -9 fe 0(ro)9 j 0(ro)9 fe 0(r o ), (5.14b) 



|V0(r o )|^ 

etc. The relationship (5.13) determines r as a function of the curvilinear coordinates £, 9 



24 



and 9?, which are orthogonal by construction. Thus 

l2 1 /<9z/,\ 2 1 (dv z \ 2 1 /<9z/ 2 \ 2 

1 /<9 2 z/ 2 \ 1 / dh e h tp \dv z - , M . W 
16 



where h%, he and h v are metrical coefficients, and <&o,<p is the operator containing deriva- 
tives only with respect to 9 and ip. It is useful to express the coefficients in (5.15) in terms 
of Cartesian coordinates: 

= |V^(r )| 2 - ^p ga i 8, ) t(r,)^(r,)9,|l(r,) + 0(( 2 ), (5.16) 

and 



d h e hp 



-A0(r) 



h^hehp \d£ h^ J 

= -A0(r o ) + O(O- (5.17) 

Inspired by our previous result (4.13) valid for the sphere, the function v z is sought in the 
limit £ — > in the form 

u z (v) = c 1 (9,i P )C + c 2 (9,ip)C 2 + --- 

+d 1 [6,ip)i ln(£/£ ) + d 2 (9, ln(£/£ ) + • • • , (5.18) 

where the formal parameter £ = i/v^ makes the argument of the logarithm dimension- 
less. The removal of diverging singularities on the rhs of Eq. (5.11), namely the vanishing 
of the coefficients attached to terms ln 2 £, ln£, £°, £ln£, £ln 2 £ and £ x in the square 
bracket, fixes 

«~wk>\- dl=0 - *=-^&- (5 - 19) 



Consequently, one arrives at 



^ ^ -wh\ ( - m^p e ™ + °« 2) ' (5 - 20) 



We present two examples. For the sphere, 0(r o ) = \Jx% + i/q + z$(= R), one has 



V0(r o ) = i^y, A0(r) = A. (5.21) 
l r o| r 



17 



Since |ro| — R and £ = R — r, one recovers (4.13) with x = £. For the cylinder infinite in 
the z-direction, 0(r o ) = \J Xq + ?/o( = one nas 

V0(r o ) = A0(r) = -L (5.22) 

l r -L| l r -L| 

where rj_ = (x ,2/o,0). Thus, 

^(r) = -e-^ 2 ln(e/6) + 0(a, (5.23) 

with £ = R- ^x 2 + y 2 . 

We would like to add that although £ has the dimension of length, it does not represent 
in general the metric distance r = |r — r |. Since one has 

T= wk)\ i+0(a (5 - 24) 

Eq. (5.20) can be expressed in terms of r as follows 

^=- t -w^)\ t2Ht/To)+0 ^ (5 ' 25) 

with the obvious definition of r . 

As was mentionned in the Introduction, v z , used through relation (1.7), is the crucial 
one-point quantity in the density-functional formalism developed in Refs. 20-22. From a 
general point of view it is desirable to have at one's disposal also off-diagonal elements of 
Green's function. The knowledge of v z is not sufficient to get the complete Green function 
within the present method (except of the case of the rectilinear hard wall). However, the 
short- distance expansion of v z (5.25), derived with the aid of curvilinear coordinates which 
mimic the global shape of the domain, depends only on the local surface characteristic. 
Since this one-point quantity is generated from Green's function itself, it is reasonable to 
expect the local shape- dependence of the latter as well. This assumption, together with 
the explicitly worked-out examples of the sphere (Sec. 4) and of the infinite cylinder, 
indicate the following boundary behavior of Green's function. For the considered pair of 
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points r, r' G f2 with the respective adjoint points ro, r' e <9f2, the perpendicular distance 
is defined by p± = |r — r' | and r — |r — r |, r' = |r' — r' |. If 7^ 0, the 3D Green function 
is supposed to exhibit the leading term oc rr' with a potential-dependent prefactor. If 
P± — 0, the analogue of the sphere result (4.9) reads 



The first term on the rhs of (5.26) corresponds to Green's function evaluated in presence 
of the rectilinear hard wall tangent to the domain surface at point r . The equation of 
the tangent plane is V0(r o ) • (r — r ) = 0. The next term is the first correction due to the 
curvature of the wall surface; the applied potential is supposed to contribute to the term 
of order O(tt'). A rigorous proof of the conjecture (5.26) is a challenge for the future. 

In conclusion, we have shown that, close to a hard wall, time- independent Green's 
functions exhibit a short- distance expansion which is to a relatively high order universal, 
i.e. does not depend on the applied potential, only the local shape of the wall is relevant. 
The universal terms are analytic in ID, and some of them are singular in 3D provided 
that the boundary curvature is nonzero. This information was obtained due to the inverse 
formalism developed in this paper, and cannot be deduced directly from the definition of 
Green's function. 
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